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Abstract 

We theoretically study quadrature and polarization squeezing in dispersive optical bistability through a vectorial Kerr cavity 
model describing a nonlinear cavity filled with an isotropic x'^"* rnedium in which self-phase and cross-phase modulation, as well 
as four-wave mixing, occur. We derive expressions for the quantum fluctuations of the output field quadratures as a function 
of which we express the spectrum of fluctuations of the output field Stokes parameters. We pay particular attention to study 
how the bifurcations affecting the non-null linearly polarized output mode squeezes the orthogonally polarized vacuum mode, 
and show how this produces polarization squeezing. 

PACS numbers: 42.50.-p, 42.50.Lc, 42.65.Sf 



I. INTRODUCTION 

Squeezing has received continued attention since the 
early eighties [2, [Sj] because it is a beautiful quan- 
tum phenomenon with amazing potential applications 
Q , the most recent of which are connected with contin- 
uous variable quantum information [^]. One particular 
aspect of squeezing that is receiving more attention at 
present is polarization squeezing, probably because its 
detection does not need homodyning, as it is the case 
with quadrature squeezing, i.e., photodetectors suffice for 
detecting polarization squeezing 0, H, 0|- Here we are 
concerned with the study of both quadrature and polar- 
ization squeezing in the field exiting a particular nonlin- 
ear optical cavity. 

One system that attracted early attention for quadra- 
ture squeezing is dispersive optical bistability [1, |^, 
[Til . Il2j | which consists in an optical cavity filled with 
a nonlinear x'^' medium which is fed with an external 
(pumping) field. For appropriate values of the system 
parameters, the output field exhibits bistability and at 
the turning points of the hysteresis cycle, perfect squeez- 
ing is obtained. Here we revisit this system because we 
are interested not only in quadrature but also in polar- 
ization squeezing. 

The generalization of this simple dispersive optical 
bistability model we consider in this article consists in 
taking into account the existence, inside the optical cav- 
ity, of a mode with orthogonal polarization with respect 
to that of the pumping field, which is also close to cav- 
ity resonance. For low pump, this extra mode remains 
empty, i.e., the output cavity field has the same polar- 
ization as the pumping field, but for large enough pump 
this orthogonally polarized mode can switch on (what is 
known as polarization instability) and then the field ex- 
iting the cavity changes from linearly to elliptically po- 
larized. This system is known as vectorial Kerr cavity 
and its mathematical modeling consists in taking 
into account self- and cross-phase modulation, as well 
as four-wave mixing, for the two orthogonally polarized 



intracavity modes, the weight of each of the three non- 
linear processes depending on the particular nonlinear 
medium under consideration 15]. 

The vectorial Kerr cavity has been studied in the past 
from the viewpoints of stability theory and pattern 
formation fisl. TTtI. [l^ [l9l. [20j . Also quantum fluctuations 
have been studied in this system but, as far as we know, 
only from the perspective of pattern formation [IJ], i.e., 
by considering a continuum of transverse modes. Here 
we consider the single transverse mode model and study 
quadrature and polarization squeezing in the field exit- 
ing the nonlinear cavity. As we show below, the vectorial 
Kerr cavity exhibits polarization squeezing. In partic- 
ular, and most interestingly, this polarization squeezing 
is optimum at the bistability bifurcation points, when 
the orthogonally polarized field is off. Thus the vecto- 
rial Kerr cavity model, which can be taken as a more 
accurate model of dispersive optical bistability, reveals 
that this well known process can act as a polarization 
squeezer. 

The rest of the article is structured as follows. In 
Section II we derive the Langevin equations for quan- 
tum fluctuations from the quantum Hamiltonian of a 
Kerr cavity pumped by a linearly polarized coherent field, 
which are derived from the Fokker-Planck equation ver- 
ified by the generalized-P distribution. In Sec. Ill, in or- 
der to keep in mind the scenario of the different solutions 
exhibited by the system, their classical stationary solu- 
tions are given and their stability properties are briefly 
reviewed. The linearized Langevin equations around the 
semiclassical stationary solutions and the spectrum of 
quantum fluctuations are introduced in Sec. IV. Sec. V 
is devoted to the analysis of quadrature squeezing. Then, 
in Section VI quantum Stokes operators and the defini- 
tion of polarization squeezing are introduced. In Sec. VII 
we state specifically the relation between the Stokes pa- 
rameters variances and the quadrature squeezing spectra 
and analyze polarization squeezing in this system, pay- 
ing special attention to the squeezing occurring at the 
bifurcations affecting the linearly polarized solution, but 



considering also polarization squeezing in the elliptically 
polarized solution in a particular case. Finally, in Sec. 
VIII we show the most relevant conclusions of our work. 



II. MODEL 

Consider an optical cavity filled with an isotropic Kerr 
medium and assume that cavity losses occur only at 
one of the cavity mirrors being the same for all possi- 
ble field polarizations (i.e., assume that the cavity is also 
isotropic). An external laser field, whose frequency is 
far enough from any resonance in the nonlinear medium 
(what makes possible that the interaction be described by 
a real cubic nonlinear susceptibility) is injected into the 
cavity. We shall assume that the injected field is linearly 
polarized, say along the x-direction. As cavity losses do 
not favour the injected field polarization, under certain 
circumstances jl6| the intracavity field can experience a 
polarization instability leading to the spontaneous gener- 
ation of a field orthogonally polarized with respect to the 
injected field (j/-polarized in our case), hence the name 
vectorial Kerr cavity for such a nonlinear system. 

The quantum Hamiltonian which describes this system 
is given, in the interaction picture, by 
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where Hfree corresponds to the free Hamiltonian of the 
two orthogonally polarized intracavity modes, H^xt de- 
scribes the injected coherent field, and Hint is the interac- 
tion Hamiltonian. In the above expressions a| and Oi are, 
respectively, the creation and annihilation operators cor- 
responding to the x-polarized (i = 1) and the y-polarized 
{i = 2) fields; loq and coc are the injected field frequency 
and the frequency of the cavity mode closest to ujq (which 
is assumed to be the same for the two orthogonally po- 
larized modes as the cavity is isotropic) , respectively; Eg 
is proportional to the amplitude of the injected coherent 
field (which we take as a real quantity without loss of 
generality); A and B are the Maker-Terhune coefficients 
(which, for isotropic media, satisfy A+B/2 — 1 15]) gov- 
erning the relative strength of cross-phase modulation 
and four-wave mixing, respectively; g — SsqHuj'^x/ ^'^^ is 
the radiation-matter coupling constant {V is the cavity 

volume, e is the medium dielectric constant and X = X 
i — x,y,z, is the nonlinear susceptibility [iBl); and, fi- 
nally, T] = ±1 takes account of the self- focusing {rj = +1 
or self-defocusing (77 = —1) cases. 
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The intracavity field exits the cavity through the sin- 
gle output mirror. Treating, as it is usual, the external 
vacuum modes as a reservoir p^ . the master equation 
governing the evolution of density matrix p of the intra- 
cavity modes is 



Ap, 
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the Liouvillian term Ap modeling the coupling between 
the system and the external reservoir through the output 
mirror (7 denotes cavity losses; 2j is the photon number 
loss rate). 



A. Fokker— Planck and Langevin equations 

In order to obtain Langevin equations describing c- 
numbers evolution from the master equation (|3ap one 
needs an appropriate quantum quasiprobability distribu- 
tion. We choose the generalized P-representation [23j . 
which is the most convenient in our case as the out- 
put field correlations are most easily derived from the 
intracavity ones with this representation: On one hand, 
it is suitable for obtaining correlations of normally or- 
dered products of operators -which provides the calcula- 
tion of fluctuations spectra outside the resonator- while, 
on the other hand, avoids the problems associated with 
other quasiprobablity distributions, such as the Glauber- 
P representation, which although also gives normally 
ordered correlation products, can exhibit non positive- 
semidefinite diffusion matrix in its dynamical equation. 

The generalized-P representation sets a correspon- 
dence between quantum operators ai{t) and aj(t) and 
independent c- numbers ai{t) and af{t), respectively, sat- 
isfying in their mean value that {ai{t))* = (^af{t)y 
We have derived the equation of evolution for the 
generalized-P distribution P (a, t) by using standard 
techniques [l3| , and have obtained the following well be- 
haved Fokker-Planck equation 
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where the new detuning parameter A is 
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with ujc = i^c — (1 + •A./2) the shifted cavity mode fre- 
quency (the shift appears because of the vacuum fluctua- 
tions contribution, which depends on the ordering chosen 
when writing the Hamiltonian [2^). Finally, the diffusion 
matrix elements are 
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being null the rest of elements. 

The Fokker-Planck equation (|4]) can be transformed 
into an equivalent set of classical-looking stochastic dif- 
ferential Langevin equation via the Ito rules [25I]. These 
Langevin equations read 



(7a 

- = A(a)+B(a).C(i): 
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where A = col (Ai, ^2, A3, A4), | = col (^,^^,^2,$^), 
with ^j(t) and ^^(t) independent white Gaussian noises 
of zero average and correlations verifying 



and matrix B (a) satisfying 
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terms are ignored). They read 
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The above equations exhibit the symmetry 
{«!, q;2, A, 77} ^ {al,a2,—A.,—ri} (remember that 
Eq has been taken to be real without loss of generality). 
Then passing from positive rj (self-focusing) to negative 
77 (self-defocusing) is equivalent to changing A ^ — A 
for what concerns the steady states and their stability 
properties. 

By introducing the following changes 
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Eqs. (|TT|) transform into an equivalent set of equations 
coinciding with those given in [1611 . Then we can use the 
results previously derived in [la| by taking into account 
the above changes. (We notice that along the following 
sections it will be more convenient to give some of the 
expressions in terms of a^, E, and r than in terms of the 
original quantities, so we shall refer later to the above 
transformations) . 

Eqs. pT|) have two different steady state solutions (see 
[iB l and references therein): The singlemode (or linearly 
polarized) solution, which corresponds to the pure Kerr 
solution [13] 
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where aj — ^JTj^e^'^j= at steady state (j = 1,2); and the 
bimode (or elliptically polarized) solution 
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III. CLASSICAL MODEL 
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Let us write the model equations ([8]) in the classical 
limit (i.e., when af is interpreted as a* and the noise 
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Notice that the bimode solution with the same intensity 
and phase <^2 = </'2s +^ ^-l^o exists, i.e., Eqs. pT|) exhibit 
phase bistability for the cross-polarized mode. 

The stability of these two steady state solutions has 
been extensively studied in il6| and turns out to be quite 
involved. Fortunately for our present purposes it will suf- 
fice to consider the stability of the singlemode solution 
and consider parameter values where the bimode solution 
does not experience secondary bifurcations, an informa- 
tion we can extract from [l6!|. From now on we shall re- 
strict ourselves to the case of liquids, for which A = 1/4 
and B = 3/2, and we refer the interested reader to [lg| 
for full details and general expressions (including possible 
anisotropics in the cavity losses). It is interesting to say 
here that there is not a threshold value for the Maker- 
Terhune coefficients for the polarization instability exists, 
what occurs is that the pump needed for reaching that 
instability tends to infinity as S ^ (a limit that corre- 
sponds to homogeneous media in which electrostriction 
is the physical origin of the nonlinearity p^). 

The singlemode (linearly polarized) solution 
(ai 7^ 0, a2 = 0), which exists for any parameter 
set, can undergo two different bifurcations: A tangent 
bifurcation (bistability bifurcation) associated with 
the appearance of optical bistability, and a pitchfork 
bifurcation (^polarization bifurcation) that leads to the 
switching-on of the orthogonal polarized field a2, i.e., 
to the bimode solution (ai ^ 0, Q!2 7^ 0). There are no 
Hopf bifurcations affecting the singlemode solution. 

The bistability of the linearly polarized solution re- 
quires that A > \/3 and the multisolution region exists 
for pump values between Eq _ and Eq _ 
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or, equivalently, for intracavity intensity values between 
hs,- and /is,-i- with 
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The polarization bifurcation, corresponding to the 
switching-on of the 02 field occurs at the pump inten- 
sity value 
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or, equivalently, at the intracavity intensity value 
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-A + V9A^ 
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For pump values larger than Eq (or lis > hs.poi) the 
solution is always the elliptically polarized one (this is 
true for liquids but not a general result [l^) and exists 
for any value of A. In Fig. 1 the two tangent bifurcations 
and the polarization instability are represented on the 
{A,E'^) plane, see Eqs.H^ 



As stated, the stability properties of the elliptically po- 
larized solution are much more involved than those of the 
singlemode solution: It can undergo secondary tangent 
bifurcations as well as Hopf bifurcations and no analyt- 
ical expressions exist for the parameter values at these 
secondary bifurcations. However these scenarios are out 
of the scope of the present work as we will consider pa- 
rameter values for which these secondary bifurcations do 
not occur. 



IV. LINEARIZED QUADRATURE SQUEEZING 
SPECTRA 

We return now to the quantum description of the sys- 
tem. In order to describe in a simple way the dynamics 
of quantum fiuctuations, the Langevin equations ([S]) are 
linearized around the stationary classical mean values of 
the intracavity fields, given by Eqs. p3)) for the single- 
mode solution (only the field parallel to the injected one 
is present in the cavity) and Eqs. ([T7)l for the bimode 
solution (the orthogonal polarized fields switches on and 
the output field is elliptically polarized). We thus write 
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with j ~ 1,2 and where the Ijs and 0^^, are, respec- 
tively, the intensity and the phase values corresponding 
to the analyzed stationary solution. Then we obtain the 
linearized Langevin equations describing the evolution 
of quantum fluctuations, 5sl = co\(Sai,5af ,5a2,Sa2), 
which read 



—Sa = A - Sa + B ■ ^(t) . 



where matrix A corresponds to 
, dA. 
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and B — B (a — as) ■ 

Quantum fluctuations are well characterized by their 
spectrum [2^, [S^ (i.e., by the Fourier transform of the 
two-time correlation functions), which are given by the 
spectral matrix M (uj) defined as 



M(w) 



dt e'"^* {5a, (0) 5a, (t)) 



(22) 



Chaturvedi et al. [2y| showed that M (u) can be directly 
obtained from Eq. ([^ 



M (uj) =(A + iujl) ^ D{A^ - iujl) ' 



(23) 



with D — BB^ and I the identity matrix. We do not give 
the explicit expression of the matrix elements of M (uj) 
because they are very lengthy. 
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Now, as we are interested in the quantum fluctuations 
of the fields outside the cavity, we must use the input- 
output theory for calculating the spectrum of quantum 
fluctuations of the output fields. As we have used the 
generalized-P representation and the cavity ispumped 
by a coherent field, the input-output theory [27| shows 
that 

/ + 00 
dt e"^' ((5a°"* (0) ^a""* (*)) = (cj) • 
-OO 

(24) 

We write down now expressions for the quadratures' 
squeezing spectra. We define the quadratures corre- 
sponding to the quantum fluctuations as 



(25) 

where we remind that j = 1,2 denotes the parallel (j = 
1) or orthogonal (j — 2) field and (5 is the (arbitrary) 
quadrature angle. 

Now we define the output quadrature squeezing spec- 
tra as 

/OO 
(:(5X°3*(t)<5X;:;^*(t-hr) :)e--dr, 
-OO 

(26) 

where : : denotes normal and time ordering. By making 
use of the input-output formalism and using the spectral 
matrix ()24p the output quadrature squeezing spectra can 
be written as 

: gr*(/3: w) : = 27 (Mffg-^^^* + M^^^e'^f^) 

-f 27 (M°2"* + M°f) , (27a) 
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It will be useful for later purposes to introduce cross- 
squeezing spectra, through the generalization of ((26)) 
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where 



ajk^Pj+Pk, Sjk^l3^-I3k. (30) 

Given the normal ordering, complete absence of fluctu- 
ations in a field quadrature means : q°^*{(},uj) := —1. 
Moreover, by taking into account the shot noise level 
outside the cavity it turns out that the symmetrically 
ordered squeezing spectra are 

g°«*(/3,c^) = 1+ :q°«*(/3,c^):, (31a) 

q",fiP„P,,iu)^ :92r(/92:/3i:^) : • (31b) 



V. ANALYSIS OF QUADRATURE SQUEEZING 
AT THE SINGLEMODE SOLUTION BIFURCA- 
TIONS 

We analyze here quadrature squeezing at the bifurca- 
tions affecting the linearly polarized solution. Following 
the steps we have outlined, after some algebra one finds 
that the quantities of interest at the bifurcations (where 
I2S — 0) can be written as 
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is twice the phase difference between the analyzed 
quadrature and the steady state phase, and jai^l is the 
normalized intensity given in Eqs. (jl2p and (|13p . In the 
above equations the frequency u has been normalized to 
7, but we do not change the symbol for not complicating 
unnecessarily the notation. 

Of course a well known general result is that at the bi- 
furcations there is a field quadrature (of the field directly 
affected by the bifurcation) that is perfectly squeezed. 
A different question is how does the bifurcation affect- 
ing one mode influence the quadrature squeezing of the 
other (orthogonal) mode. Let us consider the different 
bifurcations separately. 



A. Polarization bifurcation 

At the polarization bifurcation lis ~ Iis.poi, Eq- (|18p . 
and there is a quadrature of field Q!2 (the field that 
switches on at this bifurcation) that is perfectly squeezed 
at frequency a; = as expected (see the inset in Fig. 2). 
It can be shown that the squeezed quadrature is SX^'^p {t) 



with P = {'ippoi,opt/'2 + </'iJ and 
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which tends to zero for large negative A and to — 7r/2 for 
large positive A (Fig. 2). It is interesting to note that 
the range of -0 values for which squeezing occurs is very 
narrow. 

As for the field ai (the mode that is on) it undergoes 
large levels of quadrature noise reduction at this bifur- 
cation, even if the analyzed bifurcation does not affect 
the field ai but the field a2- These levels tend to per- 
fect squeezing for increasing positive A. In Fig. 3 the 
squeezing spectrum is shown for two values of the cav- 
ity detuning. The optimum squeezing is reached when 
■0 = 0, i.e. for (3 — (pi^, i.e., it is an amplitude squeezing, 
and occurs at the frequency 



1. Upper hraneh 

After substituting lis — Iis,+ in Eq. (j32bp one can 
show that optimum squeezing occurs for ip = n at a, fre- 
quency that is very close to zero for A close to -v/S (re- 
member that the bistability bifurcation requires A > -s/S) 
and becomes zero for A > 1.89, see Fig. 5. The squeez- 
ing level increases with detuning, as can be seen in Fig. 
6, passing from : Qj"* (""i^opt) ■— —0.75 at A = \/3 to 
: q2^'^ {-K^bJopt) '■= —0.98 for A tending to infinity. It 
is important to stress for later purposes that optimum 
squeezing occurs at = vr, i.e., for (3 = (fy^g + 7r/2. 



uJopt = y5-^A(-3A+ V8 + 9A2y. (35) 

In Fig. 4 both the optimum squeezing level and the fre- 
quency at which this squeezing occurs are represented as 
a function of cavity detuning. At this stage it is con- 
venient to remember that as A increases, the polariza- 
tion and one of the bistability bifurcations approach, see 
Fig. 1, what undoubtedly is connected with these large 
squeezing levels. 



B. Bistability bifurcations 



2. Lower branch 

After substituting lis — lis,- Eq. (j32bp one 
can show that optimum squeezing occurs for pj = tt 

at ujopt = \J [7A (A -I- VA^ ~ 3) — 15] , a frequency 
that increases linearly with A from its minimum value 
ijj — 1/a/3 at a = -^3. Fig. 7 shows the squeezing spec- 
trum for two values of the cavity detuning, and Fig. 8 
shows the optimum squeezing as well as the frequency 
at which it occurs. The squeezing level is not perfect 
(the optimum is : gf"* (t, Wopt = l/'v/S) '■— —0.75 for 
A = V3), and degrades with increasing detuning. As 
in the previous case, we stress that optimum squeezing 
occurs at = TT, i.e., for (3 = -f 7r/2. 



The bistability bifurcation has been studied several 
times in the past (see, e.g., [l^l) and we do not find it nec- 
essary to repeat here these well known results, specially 
because we shall not need them when analyzing polariza- 
tion squeezing in the following section. It will suffice to 
say that perfect squeezing, i.e. : (/3, w) : = —1, occurs 
at the bifurcations at cj = for a particular quadrature of 
the field ai. The quadratures that are perfectly squeezed 
are 
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where the subscripts "up" and "down" refer to the tan- 
gent (bistability) bifurcations of the upper and lower 
branches of the stationary solution, respectively, which 
are depicted in Fig. 1. Notice that for A 00, 
V'„p,opt and ipdown,opt I"- That is, for A cx) in 
the upper (lower) bifurcation there is intensity (phase) 
squeezing. 

Let us now analyze in more detail here the quadrature 
noise reduction at these two bifurcation in the other (i.e., 
the orthogonal) field 02, which remains off at the pump 
values at which these bifurcations occur, see Fig. 1. 



VI. QUANTUM STOKES PARAMETERS 

The polarization state of light is classically described 
by using the Stokes parameters [1^. The Hermitean 
Stokes operators are defined directly from the analogy 
with the classical parameters Q 



So = ajai 



ala2, 



S2 — a\a2 + alai, 



bi — a|ai — a2a,2, 



(38a) 



S3 — i (^a\ai — a\a2 ) . (38b) 



Operator So refers to the beam intensity, whilst Si , S2 
and 5*3 describe the polarization state. The parameter 
commutes with all the others 



So, Sj 



0, J = 1,2,3 



(39) 



whereas the remaining parameters satisfy the commuta- 
tion relation of the SU(2) Lie algebra 



Sk, Si 



^^^klva Sn 



(40) 



Therefore, the simultaneous measurements of the Stokes 
parameters are impossible in general. The non-zero com- 
mutator in (|40|) implies a restriction to the variances of 
the Stokes operators, in the form of the uncertainty rela- 
tions 

2 



VkVi > 



Sn 



I ^ m ^ I ^ 



(41) 
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where 
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Si 



(42) 



is the variance of the quantum Stokes parameter Sk ■ We 
find it convenient to normahze the variances to the mean 
intensity and thus we shaU use 



Vk 



s?.)- 



Sk 



(43) 



and then the incquaUty ([4T|) reads 

2 



VkVi > 



Sq 



I =/= m ^ k. 



(44) 



We stress that the uncertainty relations depend on 
the expected value of the operators, which makes non- 
trivial the definition of polarization squeezed states 0] ■ 
In terms of field quadratures, a squeezed state is that 
which (a) is a minimum uncertainty state (MUS) and (b) 
has fluctuations in a particular field quadrature below 
those corresponding to the vacuum state (or a coherent 
state as it is synonymous). Obviously the reduction of 
fluctuations in one field quadrature occurs at the expense 
of increasing the orthogonal field quadrature fluctuations. 
In the case of polarization squeezing this is no more the 
case because of the already commented appearance of 
the expected value in inequality ([1T|) . Let us see this is 
some more detail. A coherent state verifies @ Vfc = 1 
for all the Stokes parameters but a coherent state is not, 
in general, a MUS in terms of polarization (this depends 
on its polarization state): In effect, consider a coherent 

state such that for some m it verifies ( Sm ) < (Sa), then 



obviously VkVi > 



Sm ) / (So 



for this coherent state 



as VkVi = 1. It is also obvious that there can be states 
with some Vk < I which are not polarization squeezed 

states: Think of a state with ( Sm > = for some m and 



Vfc < 1. It is clear that is not experiencing any trans- 
fer of fluctuations from one Stokes parameter to another 
one, any rearrangement of quantum fluctuations. Then 
the quantum vacuum is not a reference for defining the 
squeezed polarization state. 

These are the reasons why polarization squeezing is de- 
fined in a different way: A field state is said to be squeezed 
if i 



V, < 



Sn 



Sc 



<Vk; I ^m^k^Q. 



(45) 



for some I, i.e., a polarization squeezed state is a state 
for which the variance of one of the Stokes parameters 



lies not only below the coherent limit, but also below 
the corresponding minimum uncertainty state limit. In a 
state verifying such a condition there is a rearrangement 
of quantum fluctuations between the different Stokes pa- 
rameters and is thus a proper squeezed state. 

We shall adopt the above criterion (|45|) but we must 
mention that there has been some debate concerning the 
most appropriate criterion for polarization squeezing and 
we refer the reader to 0, 0, [7| and references therein. 



VII. STOKES PARAMETERS' FLUCTUATION 
SPECTRA 

As stated, the polarization squeezing properties of light 
are provided by the knowledge of the variances of the 
quantum Stokes operators. As we commented in the 
previous section, polarization and quadrature squeezing 
are not equivalent, and a quadrature squeezed state does 
not imply necessarily a polarized squeezed state. Nev- 
ertheless, we find it much convenient to write down the 
spectra of the variances of the quantum Stokes param- 
eters in terms of the quadrature squeezing spectra (j26p 
and (|28|) . showing up the relation between them. We 
see below that the cross-mode spectra differentiates 
the polarization squeezing from a simple combination of 
quadrature spectra. 

After lengthy but straightforward algebra one obtains 
that the variances of the quantum Stokes operators can 
be written as 



+ \//ls/2sgj2*('^ls: '/'2s7^) 



(46a) 



2s921 



(46b) 



+ ./hJZq°2T{4>is,4>2s.^), (46c) 

^3(t^) = hsqT' (02s + f , ^) + hsqT' (01s + f , ^) , 
-^/l^sqlf (02s + f,01s + f,^) 
-Vhjisq"2l' (01s + f , 02s + f , (46d) 

where lis and correspond to the steady state values 
given by Eqs. and Eqs. ([H]) for the singlemode 

and bimode solutions, respectively. As for the quadra- 
ture squeezing spectra elements appearing in the above 
equations, they have been written accordingly to the def- 
initions in Eqs. (|26p and 
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Next we analyze the behavior of 14 (w) = Vk{ijj)/ (^SqJ, 

with (^Sq^ = lis + hs- We do this first at the bifurca- 
tions we have analyzed in Section 3, and then we analyze 
Vkiuj) as a function of pump for a particular value of the 
detuning. 

A. Polarization squeezing at the bifurcations 

Eqs. (??) simplify a lot at the bifurcation points (fT5l) 
and (jl7[) as at these particular pump values the steady 
state solution is the singlemode one, Eq. p^ . Then, 
setting l2s = one immediately obtains 



V3H 



Vi{uj) 1+ : gi((/)i3,cj) 
1+ : 92 {(Pis,^) 



= 1- 



92 



(47a) 
(47b) 
(47c) 



where we have taken into account Eq. (jSlap . 

The above expressions show that in the singlemode so- 
lution the variances of the two first Stokes parameters 
correspond to the quadrature fluctuations of the emit- 
ting (parallel) mode ai at the particular quadrature an- 
gle P = (pig (i.e., the intensity fluctuations as it must 
be) which, in general, is not the quadrature that ex- 
hibits less fluctuations. Contrarily, the two last Stokes 
parameters correspond to the quadrature fluctuations of 
the non-emitting (orthogonal) mode a2 at the particular 
quadrature angles 13 — 4>ig and /? — (p^^ + ^. 

In this case, i.e., at the bifurcations, we can give the 
explicit expressions of the Vi{uj) that read 



4|ai,|2(A-|ai,h 
ViH = 1 + ^ ^, (48a) 



6|ais|' (|ai,|' + 2A) 

1 H 

C2 



= 1 



2 lai. 



ais|' + 2A) 



V2{l^), 



where 



l-3|ai,|V4|ai,|'A. 
-4(3|aiJ^-4|aiJ'A-|-A2' 



(48b) 
(48c) 

(49a) 



C2 = 2lo 



|ai,|^ + |ai,|'A-2A^) (49b) 



ais 



(|ais 



2A 



with A ~ and \ais\^ = .g/is/7, as introduced in 

Eq. (I12a|) . and where uj is normalized to 7. Now the 
above expressions must be particularized for the different 
possible bifurcations affecting the singlemode solution. 



We must now take into account that as I2S = at 
the bifurcations, at these points the mean values of the 

Stokes parameters are (<S'o) — (Si) — hs and (S'2) = 



5*3 ) = 0, and then Eq. implies that the only Stokes 



parameters that can be squeezed are V2 and V3, and for 
that they must verify 



Si 



So 



<Vj, i,j = 2,3, t^j. 



(50) 



1. Polarization bifurcation 

In this case lis = Iis,poi, Eq- (|18p and it is not dif- 
ficult to show that neither V2 nor V3 are squeezed. We 
saw in Subsection 5.1 that optimum quadrature squeez- 
ing occurs at this bifurcation for /3 = {i'poi,optl'^ + <?^is) 
with ippoi^opt tending to for A — > —00 and to —tt/2 for 
A — > 00. It is immediately seen from this that only for 
A — » —00 it could happen that V2 be squeezed, see Eqs. 
(|47)) and Eq. (p4|) , but it turns out that this never occurs 
for finite values of the detuning because, as commented 
in Subsection 5.1, the squeezing level is extremely sensi- 
tive to the "0 value. Then there is not any polarization 
squeezing at this bifurcation. 



2. Bistability bifurcation 

In this case lis — Iis,±, Eq. The fact the op- 

timum quadrature squeezing for the field 02 occurs for 
P = (j)is + tt/2 , see Subsection 5.2, makes that V3 be 
squeezed, see Eqs. (|T7|) . Then in this particular case 
quadrature squeezing implies polarization squeezing. In 
Fig. 9 we represent the optimum (minimum) value of V3 
at the two bistability bifurcations, which occurs at the 
frequencies given in Subsection 5.2. 



B. Polarization squeezing as a function of pump 

In this subsection we study the dependence of fluctua- 
tions on pump intensity for a particular value of the de- 
tuning. We have chosen A = 1 as for this value the only 
bifurcation suffered by the singlemode solution is the po- 
larization bifurcation (remember that A > V3 for bista- 
bility exists) and thus there is not polarization squeezing 
at the bifurcation (Section VII.A.l). Furthermore, from 
[1^ we know that for A = 1 there is not any secondary 
bifurcation affecting the bimode solution. 

We concentrate on the fluctuations of the Stokes pa- 
rameters of the bimode solution. In Fig. 10 the steady 
state intensity as well as the mean values of the Stokes 
parameters are represented as a function of pump, and in 
Fig. 11 the minimum value of the variances of the Stokes 
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parameters are represented together with the mean val- 
ues of the Stokes parameters. We do not find it necessary 
to indicate at which frequency the minimum value of the 
variances is reached. Notice that polarization squeez- 
ing occurs when the represented variances fall below the 
mean value of the appropriate Stokes parameter, which 
is the one that is shown in each figure, see Eq. (|45|) . 

We see that in spite of the fact that no polarization 
squeezing exists at the bifurcation, as commented in the 
previous subsection, as pump is increased parameters Si 
and 5*2 are squeezed within some pump value ranges (we 
note that V3 increases its value for pump values larger 
than those represented in the figure). The observed po- 
larization squeezing must be attributed in this case to 
the effect of the cross-mode spectra in Eqs. (??). We 
would like to note that in spite of the fact that there is 
always some parameter for which Vi < I, see Fig. 11. 
Although this is not polarization squeezing [6i], it is a 
fact that the output field exhibits less fluctuations than 
a coherent state and this could also be useful as it is a 
form of intensity squeezing. 

These results are shown just to illustrate that the vec- 
torial Kerr cavity exhibits polarization squeezing not only 
at the bifurcations and we have checked that this hap- 
pens for other parameter values. As the behavior of the 
bimode solution can be quite involved we do not find it 
necessary to go in our analysis beyond this point. 

VIII. CONCLUSIONS 

In this article we have studied quadrature and polar- 
ization squeezing in a vectorial Kerr cavity model. After 
deriving Langevin equations for the quantum fluctuations 
from a suitable Fokker-Planck equation, namely that of 
the generalized P-distribution, we have derived expres- 
sions for the quantum fluctuations of the output field 
quadratures as a function of which we have expressed 
the spectrum of fiuctuations of the output field Stokes 
parameters. These last expressions are particularly use- 
ful for understanding the connection between quadrature 
and polarization squeezing and facilitate the study of the 
latter. 



Through the analysis of the quadrature and polariza- 
tion fiuctuations spectra, we have analyzed the condi- 
tions under which squeezing occur. We have concen- 
trated mainly on the squeezing occurring at the bifur- 
cations affecting the singlemode (linearly polarized) so- 
lution, but have considered also polarization squeezing 
in the bimode (elliptically polarized) solution in a par- 
ticular case. In particular we have demonstrated that 
polarization squeezing occurs in this system, being spe- 
cially large in the upper tangent (bistability) bifurcation 
of the singlemode solution. 

Although the results we have shown have been derived 
for particular values of the Maker-Terhune coefficients 
(namely corresponding to the case of liquids), we have 
analyzed also other cases and have convinced ourselves 
that our results are qualitatively general: Except for the 
case B — (that describes isotropic media in which the 
Kerr nonlinearity is solely due to electrostiction) in which 
the field a2 is a coherent vacuum irrespective of the pump 
value (the polarization instability does not exists in this 
case), for any B there exists polarization squeez- 
ing at the tangent bifurcations (of course the amount of 
squeezing does depend on the particular B value). 

From these results we conclude that the vectorial Kerr 
cavity model (including both cross-phase modulation 
and four-wave mixing of the two orthogonally polar- 
ized intracavity modes) constitutes a particularly inter- 
esting generalization of the standard dispersive bistabil- 
ity model from the point of view of quantum fiuctuations, 
as it predicts the existence of polarization squeezing, a 
phenomenon that, obviously, cannot be described with 
the standard model, and that exists whenever four-wave 
mixing occurs (i.e., whenever B ^ 0). It is particularly 
interesting that polarization squeezing occurs in the sin- 
glemode solution, being optimum at the tangent (bista- 
bility) bifurcations, i.e. without the need of reaching the 
polarization bifurcation. 
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Figure Captions 

Fig. 1. (a) Representation on the (i?^, a) plane of the 
two tangent bifurcations (full lines) and the polarization 
bifurcation (dashed line). The lower (down) and upper 
(up) tangent bifurcations receive these names because 
they occur, respectively, in the lower and upper branches 
of the hysteresis cycle. In (b) the output intensities are 
shown as a function of pump for A = 2. The dashed 
lines in (b) indicate unstable solutions. Another example 
of stationary solutions is shown, later, in Fig. 10(a). 

Fig. 2. The optimally squeezed quadrature of the field 
a2 (orthogonally polarized with respect to the pump) 
at the polarization bifurcation is /3 = {tppoi optl^ + '/'is)- 
V'po; opt is represented as a function of detuning. In the 



inset, the squeezing spectrum of this quadrature is rep- 
resented for A = 0. 

Fig. 3. Squeezing spectrum of the field a\ (the 
mode with parallel polarization with respect to that of 
the pump) of the optimally squeezed quadrature at the 
polarization bifurcation for the two values of the cavity 
detuning indicated in the figure. 

Fig. 4. Optimum squeezing level (full line, left verti- 
cal axis) and the frequency at which this squeezing occurs 
(dashed line, right vertical axis) for the field a\ at the 
polarization bifurcation. 

Fig. 5. Squeezing spectrum of the field (the mode 
with orthogonal polarization with respect to that of the 
pump) of the optimally squeezed quadrature at the upper 
tangent bifurcation suffered by the field a\, for the two 
values of the cavity detuning indicated in the figure. 

Fig. 6. Optimum squeezing level for the field 012 at 
the upper tangent bifurcation. 

Fig. 7. As Fig. 5 but for the lower tangent bifurcation 
suffered by the field a\, for the two values of the cavity 
detuning indicated in the figure. 

Fig. 8. Optimum squeezing level (full line, left verti- 
cal axis) and the frequency at which this squeezing occurs 
(dashed line, right vertical axis) for the field a2 at the 
lower tangent bifurcation. 

Fig. 9. Normalized variance of the Sz Stokes pa- 
rameter at the upper and lower tangent bifurcations as 
a function of cavity detuning. The optimum squeezed 
value is shown. 

Fig. 10. (a) Intensity values of the two intracav- 
ity modes and (b) normalized mean value of the Stokes 
parameters for A = 1. The polarization instability is 
marked with an arrow in the horizontal axis. The full 
and dashed lines in (b) correspond to the two possible 
values of the Stokes parameters due to the two phase 
bistability of Eqs. pTjl for the ai field. 

Fig. 11. Minimum value of the normalized Stokes 
parameters' variances (full lines) and mean values of the 
Stokes parameters (dashed lines) as a function of pump 
intensity for A = 1. Parameter Si is squeezed within 
a small pump parameter domain values around — 3, 
as can be seen in (a) and (c), while 5*1 is squeezed for 
E'^ > 5, as can be seen (c). For larger pump intensity 
values, polarization squeezing dissapears. 
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